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A METHOD FOR COMPUTING EXTREMAL 
MAXIMUM-RANGE THRUST -LIMITED ROCKET TRAJECTORIES 
WITH APPLICATION TO LUNAR TRANSPORT 

By A. Gary Childs, Ernest S. Armstrong, 
and Athena T. Markos 
Langley Research Center 

SUMMARY 

A computational technique has been described for the calculation of long-range 
extremal maximum-range thrust-limited trajectories for the transfer of a rocket vehicle 
between two points on the surface of an atmosphereless planet. The technique has been 
applied to compute such trajectories for the lunar surface with a vehicle having an initial 
lunar thrust -weight ratio of 5.0 for final to initial mass ratios from 0.95 to 0.326. Gen- 
eral characteristics are given in graphical form and were observed from the results for 
all the mass ratios. The extremal trajectories and controls are given for the particular 
mass ratios, 0.65, 0.335, and 0.326. Parameters are presented which can be used (with 
a relatively simple computer program) to reconstruct completely optimal controls and 
trajectories for the mass ratios of this study with a thrust-weight ratio of 5.0. 

INTRODUCTION 

Extensive exploration of the lunar surface may require the use of rocket -propelled 
transportation vehicles. In order to design such vehicles and plan exploration missions, 
it is helpful to know the optimal relationships between fuel and range. These relation- 
ships are especially useful when long range flights are involved. 

Manci (ref. 1) has considered the minimum fuel problem for fixed range and Childs 
and Armstrong (ref. 2) have considered the maximum range problem for fixed fuel. These 
studies involved thrust -limited rocket vehicles and assumed a uniform gravitational field 
and a soft landing. A soft landing is defined as one with near -zero velocity. Both studies 
have been applied to lunar flight. The results of these studies are valid if the resultant 
range is less than about 30 miles (48.28 km). For longer ranges both the curvature of the 
lunar surface and nonuniform gravitational field must be taken into account. Needham 
(ref. 3) has considered minimum-fuel thrust -limited trajectories with fixed range between 
points of launch and soft landing on a spherical moon with a Keplerian gravitational field. 



The longest range reported is about 65 miles (104.6 km) for a rocket vehicle with an 
initial lunar thrust -weight ratio of 2.4. Computational difficulties precluded Needham 
obtaining solutions with longer ranges. Needham noted the desirability of an effective 
computational procedure to solve such problems. 

The purpose of this paper is to demonstrate a highly effective computational proce- 
dure for obtaining optimal thrust -limited rocket trajectories in a vacuum and a Keplerian 
gravity field. This procedure combines the Pontryagin maximum principle of optimal 
control theory (ref. 4) with an initial -value algorithm for the solution of two-point 
boundary-value problems developed by Armstrong (ref. 5). Other applications of the 
procedure than those given in the present paper may be found in references 5, 6, and 7. 

The dynamic equations for the present application are presented in nondimensional 
form. The trajectories are assumed to lie in a plane containing the center of gravitational 
attraction and the rocket launch position. The Pontryagin maximum principle is applied to 
obtain necessary conditions for maximal range with specified fuel. The application of the 
maximum principle leads to a two -point boundary -value problem which is solved by the 
initial-value algorithm. The procedure is exemplified by the lunar maximum-range — 
fixed-fuel flight of a rocket vehicle with an initial lunar thrust -weight ratio of 5. A 
family of extremal solutions with payload weight as parameter and the longest ranges 
about one-half the circumference of the moon is reported herein. Also, comparisons are 
made between the theoretical necessary conditions of the maximum principle for maximum 
range with fixed fuel and minimum fuel with fixed range and it is noted where the condi- 
tions coincide. The trajectories reported are shown to qualify as minimum-fuel solutions 
between their points of take-off and landing. 

SYMBOLS 

aii elements of the matrix (i,j = 1,2,3) 

J da 

B diagonal weighting matrix 

b^ elements of the matrix B (i = 1,2,3) 

c effective exhaust velocity 

j scalar quantity measuring magnitude of e 

E energy 
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e error vector 

f s switching function (see eq. (8)) 

g gravitational acceleration 

H Pontryagin pseudo-Hamiltonian (see eq. (4)) 

I3 third-order identity matrix 

i,j indices 

J performance index 

l angular momentum 

m mass 

p^ Pontryagin auxiliary variables for maximum range problem (i = 0,. . .,5) 

r radial polar coordinate 

T thrust 

t time 

tj,t2 switching times 

W Q initial weight 

Wj- final weight 

a parameter vector 

aj elements of vector a (i = 1,2,3) 

y conditioning constant 

6a parameter correction vector 
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e eccentricity 

6 angular polar coordinate 

Aj Pontryagin auxiliary variables for minimum fuel problem (i = 0,1,. . .,5) 

l± gravitational constant 

control angle 


Subscripts: 


f final 

max maximum 

o initial 

Operations: 


(') 

o 

(■•) 

£0 

8a 

H ) 


( ) dimensionalized 
( ) differentiated with respect to time 
( ) twice differentiated with respect to time 
gradient of ( ) in the direction of the vector a 

partial derivative of ( ) with respect to a^ (j = 1,2,3) 


diag( , , ) matrix with diagonal elements ( , , ) 


matrix transpose 


= identically equal 

( J; £ y, £ j;( ) intervals on the real line 
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ANALYSIS 


Dynamic Equations 

Figure 1 shows the coordinate system and force diagram for the problem. The 
equations of motion for two-dimensional motion above an atmosphereless spherical 
planet with an inverse -square -law gravitational field are 


7 = Tc os^ + F g2_iL^ 

m r 2 

0 - X sin _ 2r 0 

r m r I 



(la) 


where 

r,6 polar coordinates 

T thrust magnitude (o = T = T max j 

*P angle between thrust and radius 

m vehicle mass 

c rocket-engine effective exhaust speed 

jl gravitational constant 

Launch from rest, soft landing on the planetary surface, and maximized range give the 
following initial and terminal conditions: 

? (tf) = ? 0 
?(t f ) = 0 

0^tfj maximized > (lb) 

J(t t ) = 0 

S(t f ) = m f 
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where the subscripts o and f indicate the initial and final values of the variables. 
Equations (la) and (lb) are nondimensionalized by use of the following relations: 



c 


m = iz- 2 — m 
Tmax 



Tmax 


c* 


where the variables on the left are nondimensional, the surface gravitational acceleration 

is g = z~7;> and the maximum thrust capability of rocket engines is T max . Nondimen- 
r 2 
L o 

sional equations corresponding to equations (la) are: 


r = l_cosJ/ + ^2 - JL 


m 


•q _ T sin i p 2r 6 
rm r 


( 2 ) 


m = -T (0 = T = 1)J 

The associated nondimensional initial and terminal conditions for equations (2) are given 
by equations (lb) without the bar notation. 


Equations of Optimal Control 


The Pontryagin maximum principle (ref. 4) is employed to obtain the equations for 
the optimal thrust magnitude and direction to maneuver the rocket from launch to a soft 
landing while $(tf) is maximized. Maximizing 0(tf) is equivalent to minimizing the 
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Lagrangian performance index 


J = 



(3) 


The problem is mathematically much like the determination of Hamilton's canonical 
equations from Hamilton's principle for conservative systems. In the present case, how- 
ever, the minimization of J is an optimization which is required in addition to the satis- 
faction of the principle and there are nonconservative forces involved. This additional 
optimization is possible because the control (a force) is left free to determine. The 
resulting equations, however, are canonically similar. The maximum principle as applied 
to the present problem states that in order for the integral J of equation (3) to be mini- 
mized, it is necessary that the following conditions be satisfied. 

The controls T (0 I T ^ 1) and \p must be chosen to maximize the so-called 
"Pontryagin pseudo -Hamiltonian" (mathematically analogous to the Hamiltonian of classi- 
cal mechanics) which is given by 


„ n ■ | a2 T cos \p u. , „ 

H = -p 0 e + Pjr + p 2 ^re + m - y - - + p 3 e + p 4 


-P 5 * 


(4) 


with the auxiliary variables (mathematically analogous to the conjugate momenta of 
classical mechanics) p^ (i = 0,1,. . .,5) satisfying 


P0 


S 0 


•2 u re T P4 sin ^ 

Pl = -P 2 e - 2p 2 -j - 2p 4 -J + — -2- 


p 2 = - Pl + 2p 4 


P3 


= 0 


H't) - °) 


(5) 


p 4 = p Q - 2p 2 r6 - p 3 + 2p 4 


P5 


= H 


p 2 cos xp + P 4 (^-^)l 
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Also, when the end time tf is not specified 


H(t) = 0 


The controls as derived in reference 8 are given by 


T = 1 


T = 0 


(*s > »)' 
(*s < «) 


(t 0 stst() (6) 


(7) 


where the switching function for the rocket engine is 

p 2 cos p 4 sin \f/ 


and, for T ^ 0, 


f s = 

® m 


rm 


P5 


( 8 ) 


or 


cos 4 / = 


P2 




sin 4/ - 


P4 


r ' P2 + 


\ry 


P4 

tan 4* = rrr- 

rp 2 


(9) 


( 10 ) 


It is shown in reference 9 that trajectories occurring with f s (t) s 0 over a nonzero 
time interval in [t 0 ,t£] (singular trajectories) cannot be optimal for minimizing fuel. 

The same proof demonstrates that singular trajectories cannot be optimal for maximizing 
range. Such trajectories will therefore not be considered. The variable Pg(t) for con- 
venience may be eliminated from the problem and replaced by f g (t) as follows. Differ- 
entiating f s (t) given by equation (8) with respect to t and making use of equations (5) 
and (9) gives 
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( 11 ) 



which is independent of p 5< 

Setting H(t) to zero at t = t Q and making use of the nondimensional initial condi- 
tions gives 


T W f s(‘o) - ■ 0 

L o 


-Pp(t 0 ) jU l 

It is assumed that fs(t 0 ) = 0. If f s (to) = 0, then — = 0. 

r o 

equation (7) T(t 0 ) = 1. In either case, 


If f s (t 0 ) > 0, then from 


*s 



P2( t o)^ 


( 12 ) 


The variable p 5 (t) may thus be eliminated and replaced by the differential equation (11) 
with initial condition (12) 

Also abnormal (p Q = 0) solutions (ref. 10, pp. 633-637) are not considered and for 
convenience Pq is taken as -1. 

From equations (7) it is seen that along a maximal range trajectory, the rocket burns 
with maximal burning rate T = 1 while f s > 0, and coasts T = 0 along a Keplerian arc 
while f s < 0. With t Q interpreted as the launch time, the rocket burns with T = 1 
from t = t 0 until such a time tf as f g changes to a negative value. The rocket then 
coasts until such a time t 2 as f s again becomes positive. This process continues 
until the first time t f for which m(t) = m f . Along the burning portions of the trajec- 
tory, the thrust direction is determined by equation (9) or (10). For the solutions pre- 
sented, it was found that over a trajectory f s changed sign twice and the trajectories 
could be characterized by initial (t e [t 0 ’^l]) ^ na -i [} € (*-2>*fj) bum periods sepa- 
rated by a single coast period (t e (tf,t 2 ^ . Then tf in terms of the specified initial 
and final mass and the two switching times is 


t f = m 0 - m f + t 2 - tf 


(13) 


When the results are collected, it is seen that extremal maximum range trajectories 
satisfy the following set of simultaneous differential equations: 
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T = 1 (f s > 0) 

T = 0 (f s < 0) 

In order to solve these equations, it is necessary to solve numerically a two -point 
boundary- value problem. 
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With the following parameters given: 


r ( t o) = r o 

* (to) = 0 


»(to) = 0 


0(t o ) = 0 
“(to) = “o 

find Pi^o), P 2 ( t o)> and P 4 (t 0 ) so that at t f 

r (tf) = r o 


t(t t ) = 0 


e(t f ) = 0 


“(tf) = m f 

The method used to solve this boundary-value problem is that discussed in reference 5 
and outlined in appendix A. The general solution of the boundary-value problem given by 
equations (14) has been programed for the Control Data 6600 electronic data processing 
system with the Fortran IV language. The program represents a useful computational 
procedure for the study of the class of thrust-limited maximum-range problems consid- 
ered in this paper. A copy of the computer program is available from the authors upon 
request. A typical family of solutions for the lunar surface is presented in the next 
section. 


RESULTS 

Extremal maximal range results are presented for a rocket vehicle with an effective 
exhaust speed, initial mass, and maximal thrust capability of 

c = 9853.2 ft/sec (3003.3 m/sec) 
m Q = 132 slugs (1926 kg) 

T m ax = 3503.8 lb (15 585.7 N) 
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Taking the lunar surface acceleration due to gravity as 


g = 5.316 ft/sec2 (1.620 m/sec2) 


gives an initial thrust -weight ratio of 


Tmax _ 5 _ 1 
m 0 g m o 


or 


m 0 = 0.2 

The value of r 0 was taken as 5.702 x 10® ft (1.738 x 10® m). 

Solutions were found for a set of values of mf/m 0 (= m|/m 0 ) obtained by fixing 
m 0 at 0.2 and varying mf so that mf/m 0 ranged from 0.95 to 0.326 in steps suffi- 
ciently small to assure convergence of the iterative procedure (see appendix A) used in 
solving equations (14). Each solution was found to have initial and final maximal thrust 
periods over the time intervals [t 0 ,tj] and (t 2 ,tfj, respectively, separated by a coast 
period over (tj^tjj]. Results needed to reconstruct the trajectories are given in table I. 
Solutions were not obtained below mf/m 0 = 0.326 because of convergence difficulties in 
the iterative procedure. Past this point the terminal values of r, r, and 6 became 
sufficiently insensitive to changes in the p^(t 0 ) (i = 1,2,4) to lower the convergence rate 
to a point where the acquisition of additional solutions was too expensive in terms of com- 
puter time. This difficulty is reflected by the sharp increase in the values of the auxiliary 
variables at mf/mQ = 0.326. 

TABLE I.- CHARACTERISTICS OF EXTREMAL MAXIMUM RANGE SOLUTIONS 


Mass ratio 

Initial values of auxiliary variables, 
nondimensional 

Switching and final times, 
sec 

m f/ m 0 

pita) 

P2(‘o) 

P4 ta) 

*1 

h 

tf 

0.95 

0.01096 

0.05603 

0.01884 

9.414 

60.48 

69.62 

.89 

.03552 

.1288 

.04323 

21.08 

138.9 

158.6 

.83 

.07606 

.2096 

.07041 

33.19 

225.3 

255.3 

.75 

.1647 

.3355 

.1134 

50.07 

357.7 

400.4 

.65 

.3683 

.5433 

.1873 

72.71 

566.6 

624.0 

.55 

.7797 

.8569 

.3096 

96.92 

856.9 

927.0 

.45 

1.814 

1.455 

.5923 

123.5 

1329 

1415 

.40 

3.178 

2.072 

.9712 

138.0 

1738 

1827 

.38 

4.170 

2.448 

1.256 

143.6 

1954 

2041 

.37 

4.932 

2.706 

1.481 

146.6 

2085 

2172 

.36 

6.000 

3.034 

1.802 

149.7 

2233 

2321 

.35 

7.650 

3.478 

2.311 

152.7 

2404 

2493 

.335 

13.81 

4.758 

4.278 

157.2 

2723 

2813 

a .326 

74.94 

13.49 

24.51 

159.7 

3000 

3090 


a Negative altitude during early part of trajectory. 


Range, 
miles (km) 


r o^f) 


1.973 

( 3.175) 

10.27 

( 16.53 ) 

26.59 

( 42.79 ) 

64.96 

( 104.5 ) 

153.8 

( 247.5 ) 

321.1 

( 516.8 ) 

670.6 

(1079 ) 

1024 

(1648 ) 

1229 

(1978 ) 

1365 

(2197 ) 

1530 

(2462 ) 

1738 

(2797 ) 

2216 

(3566 ) 

2938 

(4728 ) 
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In all cases convergence was obtained to within ±I0 - 3 feet (3 x 10"^ m) of 
r(tj) = r Q) ±0.1 ft/sec (0.03 m/sec) of r(t^) = 0, and ±10"® rad/sec of 0 (tf) = 0. 

For each solution, Ps(t^ > 0 whereby it is also a minimum fuel extremal. (See 
appendix B.) 

Figure 2 gives lunar surface range r Q 0 (tf) as a function of final to initial mass 
ratio Surface range appears to increase exponentially for decreasing payload 

nif and thus indicates that for larger amounts of fuel, the extremal trajectories deliver 
many more units of range per unit mass of fuel. 

As can be seen from table I, the duration of the coast periods increases for 
decreasing payload varying from about 74 percent of the total flight time for 
m^m 0 = 0.95 to about 92 percent for nif/m 0 = 0.326. For all values of mf/m 0 , the 
range accrued over the coast period was found to be the predominant part of the total 
range and, overall, the range accrued over the coast periods increased with increasing 
fuel. The percentage of fuel used for injection into the coast phase increased with 
increasing fuel supply and ranged from 3.7 percent (nif/m 0 = 0.95) to 43 percent 
(mf/m 0 = 0.326). 

For all mf/m 0 the coast orbits were elliptic and became nearly parabolic for 
mf/m 0 = 0.95. Figure 3 shows the coast phase eccentricity as a function of m^m 0 . 

The trajectories displayed increasing coast energy and angular momentum for 
increasing fuel supply. Maximal radial distances on the trajectories occurred during 
the coast periods. Maximal altitude r max - r 0 as a function of surface range r Q 0 (tf) 
is shown in figure 4. The points of maximal radial distance are also the points of mini- 
mal 0 . Figure 5 shows that the minimal 0 at the peak altitudes on the coast orbits 
increases with fuel supply. The behavior of r max - r 0 in figure 4 may be explained 
by recalling that over the coast periods, 



and noting that, for given angular momentum T and mass in, minimal 9 may be 
increased by decreasing r max . However, the case of a low orbit is tempered by the 
fact that lower orbits possess higher kinetic energy which must be dissipated to make a 
soft landing. Thus, lowering the orbit is postponed until the payload is sufficiently small 
to give a large deceleration capacity for soft landing. 

Figure 6 presents detailed results for nif/m 0 = 0.65. Figure 6(a) gives the thrust 
angle if/ over the launch and landing phases. The angle if/ is undefined over the coast 
phase. Figure 6(b) gives altitude r - r Q as a function of range r o 0(t). Figure 6(c) 
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gives altitude as a function of flight time. Figures 6(d) and 6(e) give radial r and 
tangential r£> velocity time histories. Figure 7 gives similar results for 
m f /m 0 = 0.335. 

The general characteristics given in figure 6 for mf/m 0 are typical of the inter- 
mediate range solutions. Results for a typical long-range solution m.0n o = 0.335 are 
given in figure 7. Figure 8 presents results for mf/niQ = 0.326, the smallest ratio for 
which a solution was obtained. For this case a dip beneath the lunar surface is noted 
during the initial thrusting phase. Such dips will always occur when the initial angle 

m g 

90° - ^ is less than ll°32 t = sin“l =r -^ — and causes the radial component of thrust 

^max 

to be less than the gravitational force. 

From figure 3 it appears that for some value of mf/m 0 less than 0.326 (approxi- 
mately 0.320 by extrapolation), the coast trajectory will become circular. Figure 9 
shows coast orbit eccentricity as a function of maximal altitude. This curve is not as 
easily extrapolated as that of figure 3 to estimate the altitude of the circular orbit. 
However, the specific energy (energy per unit mass) of the coast orbits is monotoni- 
cally increasing for decreasing mf/m 0 with E/m = -1.5291 x 10 7 ft2/sec 4 
(-0.1421 x 10 7 m2/sec4) for m f /m 0 = 0.326 and E/m = - JL = -1.5154 x 10 7 ft2/sec 4 

(-0.1408 x 10 7 m2/sec 4 ) for a circular orbit about the lunar surface. It is thus con- 
jectured that the circular orbit will have near-zero altitude. Since no altitude con- 
straints have been theoretically imposed, the thrusting part of this trajectory may con- 
tain negative altitudes as did the mf/m 0 = 0.326 case. 

If it is recalled that the range accrued over the coast phase was, in all solutions 
obtained, the predominant part of the total range, it is evident that the most important 
characteristic of each extremal is the coast phase. This result should follow for other 
lunar rocket vehicles as well as for the one considered herein. Although, in practice, it 
may not be possible to follow exactly an extremal trajectory for a given vehicle and value 
of m£/m 0 , an attempt could be made to achieve the corresponding coast phase as effi- 
ciently as possible. In this way near theoretical extremal range should be obtained. 

Solutions for rocket-engine configurations different from the one considered herein 
may be similarly obtained by the computational procedure presented. 

CONCLUDING REMARKS 

A program for obtaining maximum- range thrust-limited rocket trajectories between 
points of launch and soft landing on the surface of an atmosphereless planet has been 
described. An application of the program for the lunar surface has produced trajectories 
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for a vehicle with a lunar thrust- weight ratio of 5.0. Results are presented for final to 
initial mass ratios from 0.95 to 0.326, the longest trajectory being about 3000 miles 
(4827.9 km). The trajectories presented were shown in an appendix to qualify as mini- 
mum fuel extremals between the points of launch and soft landing. 

Since the primary purpose of this paper is the demonstration of the capability of 
calculating long maximum range trajectories in a Keplerian gravitational field, results 
have only been obtained and presented for one value of lunar thrust-weight ratio. How- 
ever, two facts are clear from the results given. First, as fully expected, the range 
accrued during the coast period was found to be the predominant part of the total range. 
Thus, it may very well be possible to develop excellent suboptimal trajectories by using 
the same coasting phase with more practical thrusting phases (especially for the "dip" 
trajectories). Secondly, for larger amounts of fuel, the extremal trajectories deliver 
more units of range per unit mass of fuel; that is, for longer ranges, the trajectories are 
increasingly efficient. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 27, 1970. 
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APPENDIX A 


SOLUTION OF THE BOUNDARY- VALUE PROBLEM 


The approach taken to obtain solutions of the boundary-value problem represented 
by equations (14) is an iterative method discussed in reference 5. In reference 5, for a 
similar boundary- value problem, a vector e(<5,tf) is defined so that when e («,tf) = 0, 
the terminal boundary conditions are satisfied and the unknown initial conditions are a 
where the arrows indicate vector quantities. For the problem at hand, 



ei(a,tf) 


r(a,tf) - r Q 

e(o,tf) = 

e 2 (“>tf) 


r(a,t f ) 


^(^tf) 


0(a,tf) 


(Al) 


and 



“i 


Pl(*o) 

a = 

“2 

= 

P2( t o) 


_“3_ 


_P4( t o) 


(A2) 


The argument a is included in r(c?,tf), r(c?,tf), and 9(3 ,tf) to indicate their implicit 

dependence on a through equations (14) and also serves to indicate that a specific a 
was used to obtain the value of the variable at tf. 

The magnitude of e"(a,tf) is measured by a scalar quantity 


E[e(3 >tf )] = l 


(A3) 


where B is a positive definite diagonal matrix of weighting elements and primes 
denote a matrix transpose. Here 


E [e(a,tf)] = |jbi[r(a,tf) - r 0 ] ? + b 2 r 2 (a,tf) + 


where 


(A4) 


B = diag(b 1 ,b 2 ,b 3 ) 
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APPENDIX A - Continued 


Initially, a value of a is assumed, the differential equations (14) are integrated 
forward in time until the first time at which 


m(t) - mf = 0 


and e(a,tf) is evaluated. If e(a,tf) vanishes, or, for numerical purposes, is suffi- 
ciently small, the boundary-value problem is considered to be solved. Otherwise, the 
assumed a is corrected by 


6a = - 


9e 


K S,tf )' B H( S ’ t f) +yI 3 


9a 


i|(a,t f )’ Be(a,tf) 


(A5) 


where the superscript -1 denotes a matrix inverse, I3 is a 3 by 3 identity matrix, 
and y > 0 is adjusted so that 

EjV(a +6a,tf + 6t^|j < E|lr(a,tf)J (A6) 

The increment 6tf is added to tf in the expression (A6) since the final time for 
which m(t) - mf = 0 by using a may not be the same as that obtained by using a+ 6a. 

For this problem ^(a,tf) is a 3 by 3 matrix with elements a^ (i = 1,2,3; j = 1,2,3) 

given by 


“21 ■ ) 

3 

90 /— , \ 

a 31 " 

j J 


(A7) 


The quantities in equations (A7) are obtained by finding simultaneous solutions at 
tf of equations (14) and a set of sensitivity equations derived as shown in reference 5 
(p. 21, case (3c)) with tf determined through 


m(t) - mf = 0 


With the correction 6a computed, the assumed a is replaced by a + 5a and the 
process is repeated by treating a + 6a as the assumed value. 
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APPENDIX A - Concluded 


For the results of this paper, initial assumed values of a were obtained by using 
a small amount of fuel and the results of reference 2 for the shorter trajectories and 
then using the a from the last obtained solution as an initial guess for larger amounts 
of fuel. In each case 


B = diag(l,l,l) 

This approach for obtaining assumed values for a is also suggested for values of initial 
thrust-weight ratio other than 5.0. 
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APPENDIX B 


RELATIONSHIPS BETWEEN THE MAXIMUM RANGE 
AND MINIMUM FUEL PROBLEMS 

In this section comparisons are drawn between the necessary conditions of the 
Pontryagin maximum principle as applied to the maximum-range — fixed-fuel and 
minimum-fuel — fixed-range problems. Each problem is for a rocket vehicle launched 
from rest and executing a soft landing on the surface of a spherical planet with an 
inverse-square-law gravity field. 

The necessary conditions of both problems are summarized in the following chart. 
Conditions common to both problems are centered in the chart. The auxiliary variables 
of the maximum principle are denoted by pj for the maximum-range — fixed-fuel prob- 
lem and by for the minimum-fuel — fixed-range problem. 


NECESSARY CONDITIONS FOR THE MAXIMUM RANGE AND MINIMUM FUEL PROBLEMS 


Maximum Range-Fixed Fuel 


Minimum Fuel-Fixed Range 


The nondimens ional dynamic equations for both problems are 
- = T C ° 3 * + r § 2 - n/r£ (r(t 0 ) = r(t f ) = r Q , r(t D ) = r(t f ) = oj 

(e(t 0 ) = S ( t o) = ®(tf) = °) 

(m(to) = “o) 

where 0 < T < 1 , t f unspecified 


0 = T sin t 2r9 


m = -T 


( 0 (t f ) maximized, m(tf) specified) 


( 9 ( t-f ) specified, dIq - m^* minimized) 


In integral form minimize 
't, 


r - - 


r-* 


mdt 


Application of the Pontryagin Maximum Principle leads to the 
necessary conditions 


T = 


(fg > 0) 
(*8 < 0 ) 
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APPENDIX B — Continued 


where 
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m o 

r 


+ * e + J U*L* . Sf 

3 rm r 


- ft - \) 1 ■ ° 


Given an extremal maximum range solution with specified fuel expenditure m 0 - mf 
and P 5 (tf) ^ 0, this same solution can be made to satisfy the necessary conditions for 
fuel optimality between the maximum range launch and landing points with extremal 
fuel given by m 0 - mf by maldng the identification 


M(t) = PfW (i = 1,2,4) 


* 3 ( 1 ) = P 3 W ~ Pq 


^■5(1) “ *-Q = P5W 
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APPENDIX B - Concluded 


In order to satisfy Ag(tf) = 0, let 

X 0 = “Psfc) 

Also Ag(tf) = -pq follows from pg(tf) = 0. 

Likewise, if given an extremal minimum fuel solution with extremal final mass mf, 
specified range angle 0(tf), and ^(tf) g o, this same solution can be made to satisfy the 
necessary conditions for maximal range with specified final mass mf and extremal 
range angle 9 (tf) by making the identification 

Pi(t) = Xf(t) (i = 1,2,4) 


P 3 W - Pq = * 3 (t) 


P5(t) - ^s(t) ” ^0 


In order to satisfy pg^tf) = 0, let 


Pq = 


Each of the maximum range extremals presented as results in this paper had 
Pg(tf) > 0 whereby they are minimum fuel extremals. 

The program may thus be used to construct minimum fuel extremals by finding a 
maximum range extremal which achieves the desired range and then checking to deter- 
mine whether Pcftf)>0. 
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Figure 6.- Optimal control and trajectories. mf/m 0 = 0.65. 
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(e) Variation of tangential velocity with time. 
Figure 6.- Concluded. 
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(a) Variation of thrust angle with time. { i}j is measured from the vertical direction as shown in fig. 1.) 
Figure 8.- Optimal control and trajectories. mfy'm 0 = 0.326. 
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Figure 8.- Concluded. 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 


03U 001 55 51 3DS 00903 

AIR FORCE WEAPONS LABORATORY /WLOL/ 
K.IRTLAND AFB , NEW MEXICO 87117 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


ATT E* LOU BOWMANf CHIEF*TECH. LIBRARY 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


r T he aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knoivl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 

TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 

Details on the availability of these publications may be obtained from : 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPEQAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS : Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


